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ABSTRACT 

Aims. We observed the high energy peaked BL Lac object PKS 0548-322 (BZB J0550-3216) with Swift to study the temporal and spectral 
properties of its synchrotron emission simultaneously in the Optical, Ultraviolet and X-ray energy bands. 

Methods. We carried out a spectral analysis of 5 Swift XRT and UVOT observations of PKS 0548-322 over the period April - June 2005. 
Results. The X-ray flux of this BL Lac source was found to be approximately constant at a level of Fp-io keV) — 4 X 10~" erg cirr 2 s -1 , a 
factor of 2 brighter than when observed by BeppoSAX in 1999 and close to the maximum intensity reported in the Einstein Slew Survey. The 
very good statistics obtained in the 0.3-10 keV Swift X-ray spectrum allowed us to detect highly significant deviations from a simple power 
law spectral distribution. A log-parabolic model describes well the X-ray data and gives a best fit curvature parameter of 0.18 and a peak 
energy in the Spectral Energy Distribution of about 2 keV. The UV spectral data from Swift UVOT join well with a power law extrapolation of 
the soft X-ray data points suggesting that the same component is responsible for the observed emission in the two bands. The combination of 
synchrotron peak in the X-ray band and high intensity state confirms PKS 0548-322 as a prime target for TeV observations. X-ray monitoring 
and coordinated TeV campaigns are highly advisable. 



Key words, radiation mechanisms: non-thermal - galaxies: active - galaxies: BL Lacertae objects - X-rays: galaxies: individual: PKS 
0548-322 (BZB J0550-3216) 



1. Introduction 

BL Lacertae objects (BL Lacs) constitute a rather peculiar and 
extreme class of Active Galactic Nuclei (AGNs). The main 
characteristics that distinguish BL Lacs are the rapid variability 
at all frequencies, high and variable radio and optical polariza- 
tion, compact and flat-spectrum radio emission, superluminal 
motion, smooth and broad non-thermal continuum covering the 
electromagnetic spectrum from radio to y-rays and the almost 
complete absence of emission lines in the optical band. The ex- 
treme properties of BL Lacs are successfully explained in terms 
of relativistic beaming, i.e. of a relativistic bulk motion of the 
emitting region toward the observer dBlandford & Rees 1978 
|Urry & Padovani 1995) . 
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The Spectral Energy Distribution (SED) of BL Lacs is gen- 
erally characterized, in a Log{v F(v))-Log{v) representation, 
by two emission peaks: the first is produced by synchrotron 
emission by relativistic electrons in a jet closely aligned to the 
line of sight, while inverse Compton scattering by the same 
population of relativistic electrons is thought to be at the origin 
of the higher energy peak (e.g. Ghis ellini & Maraschi 19 89 1. 
BL Lacs are often divided into two classes according to the 
position of the synchrotron energy peak: low energy peaked 
BL Lacs (LBLs), with the peak located at IR/optical wave- 
lengths, and high energy peaked BL Lacs (HBLs) with the 
synchrotron emission peaking in the UV/X-ray energy band 
(Giommi & Padovani 1994; Padov ani & Giommi 19951 1. 

PKS 0548-322 (z=0.069, |Fo5bury & Disney 1976) , also 
named BZB J0550-3216 in the recent Multifrequency 
Catalogue of Blazars (Mas saro et al. 20051 . is a remarkable 
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Table 1. Swift observations and exposures of PKS 
0548-322 in spring 2005. 



Date 


Start UT 


XRT Exp. 


UVOT Exp. 






(s) 


(s) 


April 1 


00:43 


1,258 


/ 


April 26 


23:31 


5,212 


7,791 


April 28 


10:53 


1,354 


/ 


May 13 


14:22 


3,362 


3,380 


May 21 


10:09 


9,224 


9,135 


May 22 


22:52 


40,191 


/ 


May 24 


00:02 


1,349 


5,402 


May 26 


00:46 


400 


1,926 


May 29 


00:07 


8,083 


8,715 


June 14 


16:23 


1,235 


/ 


June 24 


17:15 


7,951 


/ 



BL Lac object characterized by a relatively strong and fast 
variability in the X-ray energy band (Bl ustin et al. 2004b . It 
is hosted in a giant elliptical galaxy (Falomo et al. 1995, 
IWurtz et al. 19 96 1 which is the dominant member of a rich 
cluster of galaxies. The synchrotron power peaks in the X- 
ray band and for this reason it is classified as an HBL 
source JPadovani & Giommi 19951 . PKS 0548-322 was ob- 
served on several occasions by a number of X-ray astronom- 
ical satellites showing strong variations both in intensity and 
in spectral shape. A brief summary of historical X-ray data 
on PKS 0548-322 can be found in |Costamante et al. (2001)| 
who also presented the analysis of three BeppoSAX obser- 
vations. The brightest flux was reported in the Einstein Slew 
Survey (Perl man et al. 19 96 ). The possibility that the X-ray 
spectrum of PKS 0548-322 might deviate from a sim- 
ple power law was apparent since the very early observa- 
tions. |Urry et al. (1986)| using Einstein data, interpreted the ob- 
served curvature as due to an excess low energy absorption, 
whereas Madejski (1985) suggested a spectral bending. PKS 
0548-322 was observed by EXOSAT at five epochs between 
1983 and 1986 JBarretal. 19881 |Garilli & Maccagni 1990| l 
and broken power law best fits gave a peak energy of the 
Spectral Energy Distribution in the 2.5-5 keV range. Similar 
results were obtained in the re-analysis of all EXOSAT obser- 
vations by |Ghosh & Soundararajaperumal (1995)| who noticed 
that a log-parabolic law represents well the multifrequency 
spectrum of PKS 0548-322. A subsequent GINGA observa- 
tion in February 1991 (Ta shiro et al. 19951 showed a much 
flatter spectrum, with an X-ray photon index of 1.84+0.02 
in the 2-30 keV range implying that the synchrotron peak 
in the SED had moved at energies larger than 30 keV. PKS 
0548-322 was observed twice by BeppoSAX in 1999 on 
February 20 and on April 7. A first analysis was reported 
by |Costamant e et al. (200 1)| who found evidence of intrinsic 
curvature because a broken power law model was necessary 
even in presence of an extra absorption. The possibility of 
a presence of circumnuclear ionized gas was suggested by 
|Sambruna & Mushotzky (1998)| who reported evidence of an 



- 


0.3-1.5 keV 

—\ 1 1 1 


■ 1 1 1 1 1 1 1 1 — 


.5-10 keV ; 


: 1 1 1 1 1 1 1 1 — 

- HR = (1.5 


— | 1 1 1 

-10)/(0.3-1.5) Z 



5x10 4 10 5 

Time (s) 



Fig. 1. Swift XRT 0.3-1.5 keV (upper panel) and 1.5-10 keV 
(middle panel) light curves of PKS 0548-322. In the lower 
panel the corresponding Hardness Ratio (HR) is plotted. Data 
are binned to 12 ks resolution and error bars indicate statistical 
uncertainties at the lcr level. 

absorption feature around 0.6 ke V with a neutral hydrogen col- 
umn density Nh — 10 21 cm -2 , however this feature was not 
confirmed by a subsequent spectroscopic XMM-Newton obser- 
vation (Blusti n et al. 200 4 1. Finally, the Chandra detection of a 
diffuse soft X-ray radiation around PKS 0548-322, interpreted 
as a thermal emission of the host galaxy on kpc scale, was re- 
ported by Dona to et al. (2003)1 

This interesting source was observed on several occasions 
by the Swift satellite ( Gehr els et al. 20 04 1 from April to June 
2005. In this paper we present the results of a detailed spectral 
analysis of the X-Ray Telescope (XRT, IBurrows et al. 2 005 1 
and Ultraviolet/Optical Telescope (UVOT, |Roming et al. 2005) 
data, confirming that the X-ray spectrum of this HBL object 
shows a well established curvature similar to that found in other 
sources of the same type. In Section |2 the observations and 
the data reduction are presented, in Section [3] we describe the 
Swift XRT spectral analysis and Section |4] is dedicated to the 
BeppoSAX data analysis. Finally the results are discussed in 
Section |5] Throughout this paper errors are quoted at the 90% 
confidence level for one parameter of interest (A^ 2 = 2.7) un- 
less otherwise specified. 

2. Observations 

As part of a Swift key project dedicated to the observation 
and monitoring of a sample of Blazars PKS 0548-322 was 
pointed eleven times over the period April-June 2005. The 
journal of these observations is given in Table ^ where we 
also report the net exposures with the XRT and UVOT instru- 
ments. On May 22 a deep (about 40 ks) XRT exposure of PKS 
0548-322 was taken. The exposure times in the Burst Alert 
Telescope (BAT, B arthelmy et al. 2005) were not sufficient to 
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Table 2. Best fit spectral parameters of the log-parabolic model for the XRT (first section) and BeppoSAX (second section) 
observations of PKS 0548-322. Numbers in parenthesis are statistical errors at the 90% confidence level. Fluxes in the 2-10 keV 
energy band are corrected for Galactic absorption. 



Instrument 


Date 


K 


a 


b 


E P 


■^2-10 


A^/dof 












(keV) 


(erg cirT 2 s~') 




Swift XRT 


2005 April 26 


1.41 (0.05) 10~ 2 


1.81 (0.07) 


0.18(0.13) 


3.4 


4.0x10-" 


0.75/118 


Swift XRT 


2005 May 21 


1.48 (0.04) 10" 2 


1.71 (0.05) 


0.38 (0.10) 


2.4 


4.0x10-' 1 


0.90/196 






1.46(0.04) 10~ 2 


1.78 (0.03) 


0.18 frozen 


3.5 


4.2xl0~ n 


0.94/197 


Swift XRT 


2005 May 22 


1.66 (0.02) 10- 2 


1.90 (0.02) 


0.18 (0.04) 


1.9 


4.1x10-" 


1.03/408 


Swift XRT 


2005 May 29 


1.31 (0.04) 10~ 2 


1.77 (0.06) 


0.34 (0.11) 


2.2 


3.3x10"" 


1.01/181 






1.29 (0.04) 10~ 2 


1.83 (0.04) 


0.18 frozen 


2.7 


3.5x10-" 


1.03/182 


Swift XRT 


2005 June 24 


1.52 (0.05) 10~ 2 


1.79 (0.06) 


0.26 (0.11) 


2.5 


4.1x10-" 


1.09/183 


BeppoSAX LECS+MECS 


1999 Febr. 20 


0.77 (0.07) 10~ 2 


1.57 (0.10) 


0.48 (0.10) 


2.8 


2.3x10-" 


1.23/55 


BeppoSAX LECS+MECS 


1999 April 07 


0.69 (0.06) 10 2 


1.82 (0.10) 


0.42 (0.10) 


1.6 


1.5x10-" 


0.91/55 



detect a source with a typical intensity lower than 2 mCrab like 
PKS 0548-322. 



2.1. XRT data reduction 

To obtain a signal with sufficient statistics for a detailed spectral 
analysis we only considered observations longer than 5 ks. In 
particular, the X-ray spectrum accumulated during the longest 
XRT observation of May 22 has an excellent photon statistics 
and allowed us to perform a very accurate spectral analysis. 
All XRT observations were carried out using the most sensi- 
tive Photon Counting readout mode (see IHill et al~ 004 for a 
description of readout modes). The XRT data were processed 
with the XRTDAS software package (v. 1 .8.0). Event files were 
calibrated and cleaned with standard filtering criteria with the 
xrtpipeline task using the latest calibration files available in the 
Swift CALDB distributed by HEASARC. Events in the energy 
range 0.3-10 keV with grades 0-12 were used in the analysis 
(see lBurrows et al. 2005l for a definition of XRT event grades). 

The source count rate was high enough to cause some pho- 
ton pile-up in the inner 6 pixel (~ 14") radius circle within 
the peak of the telescope Point Spread Function (PSF), as de- 
rived from the comparison of the observed PSF profile with 
the analytical model reported by Mor etti et al. (2005)| We thus 
avoided pile-up effects selecting events within an annular re- 
gion with an inner radius of 6 pixels and an outer radius of 30 
pixels. The background was extracted from a nearby source- 
free circular region of 50 pixel radius. Ancillary response files 
for the spectral analysis were generated with the xrtmkarf task 
applying corrections for the PSF losses and CCD defects. The 
latest response matrices (v. 008) available in the Swift CALDB 
were used. The spectrum was binned to ensure a minimum of 
20 counts per bin, and energy channels between 0.4 keV and 
0.6 keV were excluded to avoid undesired effects on the mea- 
sured spectral parameters due to residual instrumental features 
(Campana et al. 2006). 



2.2. UVOT data reduction 

UVOT observations with good exposure times were available 
for a number of pointings (see Tableful. Sky corrected images 
were derived from the Swift archive and aperture photometry 
was made with UVOTSOURCE using a 6" (12 pixels) radius 
for the V,B, U filters and 12" for the W\,M2 and W2 filters. 
Count rates in all the filters were always less than 5 counts/s, 
well below the pile-up threshold (10-15 counts/s). 

The source showed no appreciable variations in all the 
pointings, with an average value < B >= 16.9. We report here 
only the data regarding the May 21 observation, which is the 
longest and is nearly simultaneous with the longest XRT ob- 
servation of May 22. 

The host galaxy of this source is a giant elliptical that has 
been extensively studied (e.g. by Falomo et al. 1995). To es- 
timate the galaxy contribution within our 6" aperture we in- 
tegrated the best fit De Vaucoulers profile published by them 
(their Fig.4) obtaining R-15.7. Assuming the typical colors 
for elliptical galaxy (B - V=0.96, V - fl=0.61; Fukugita et 
al. 1995) we derive V=16.3 and B=17.2. The AGN luminosity 
is therefore just 25% of the observed flux in these bands. For 
this reason we considered only photometric data in the U and 
UV bands where the host galaxy contribution is smaller. The 
measured magnitudes in the U and UV filters were U = 16.5, 
UVWl = 16.4, UVM2 = 16.7, UVW2 = 16.6. As the zero 
point uncertainty in the optical bands of UVOT is about 0.1 
mag (about 10% in flux) an accurate estimate of the spectral 
shape is not currently feasible. 

3. The X-ray spectrum 

To study the X-ray spectral distributions of PKS 0548-322 we 
first considered the deep 40 ks exposure performed on May 22. 
In Fig.^ we show the 0.3-1.5 keV (upper panel) and 1.5-10 
keV (middle panel) light curves of PKS 0548-322, together 
with the corresponding hardness ratio (lower panel). From the 
figure is apparent that no significant temporal and spectral vari- 
ability was present during the observation. We thus performed 
the spectral analysis using the events accumulated in the en- 
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tire duration of the observation. We adopted the following two 
spectral models, a single power law: 

F(E) = K E~ a photons cirrVkeV" 1 (1) 

and a log-parabolic law ( Mass aro et al. 2 004a i: 

F(E) = K E- {a+bLogE) photons cirrV'keV" 1 (2) 

where the parameter a is the photon index at 1 keV while b 
measures the curvature of the parabola. The latter model has 
the property of describing curved spectra with only one addi- 
tional parameter with respect to the single power law. The log- 
parabolic model is also very useful for the estimation of the 
energy and the flux of the SED peak, simply given by, respec- 
tively: 

E p = \0 (2 - a)/2b keV, (3) 
v p F(v p ) = (1.60xlO- 9 )K10 {2 - a)2/4b ergs" 1 cm" 2 . (4) 

The log-parabolic law generally fits well the X-ray spectra of 
HBL sources, as shown by Massaro et al. (2004a,b) for the 
cases of the X-ray spectra of Mkn 421 and Mkn 501. The 
relations between the parameters of synchrotron and inverse 
Compton radiation emitted by an electron population with a 
similar energy distribution have been recently investigated by 
Massaro et al. (2006). 

When considering curved spectra a parameter that needs 
special attention is the amount of low energy neutral hydrogen 
absorbing column density A^h- We first fitted the background 
subtracted 0.3-10 keV spectrum with a single power law model 
with the absorption column density fixed at the known Galactic 
value of jV H =2.49xl0 20 cm" 2 (|Murphy e t~l. 1996| >. We ob- 
tained a photon index of 1.96+0.02 with reduced;^ 2 of 1.15 
with 409 degree of freedom (dof). From the inspection of 
the residuals a clear systematic effect reflecting the presence 
of spectral curvature was observed. We thus adopted the log- 
parabolic model with the same amount of A^h- The model pro- 
vided a very good fit with;^ 2 =1.03 (408 dof), and a;^ 2 decrease 
of 50 for only one additional parameter. The F test gives a prob- 
ability of about 10" 11 that this improvement is due to chance. 
Figure [2] shows the best fit spectrum with this model and the 
residuals. 

The same log-parabolic spectral law was then also applied 
to the analysis of the observations of shorter duration. The re- 
sulting best fit parameters, the SED peak energies E p and fluxes 
in the 2-10 keV band are given in Table |2 We note that for 
these shorter observations the parameters a and b are charac- 
terized by statistical uncertainties much larger than those of the 
May 22 observation. In two cases (April 26 and June 24) the 
value of b was found to be consistent with the one measured on 
the May 22 deep observation, while on May 21 and 29 it was 
larger (b ~ 0.3 - 0.4). For these two latter pointings, best fits 
with the b value frozen at 0.18 (the value May 22 value) gave 
very small increases of the^ 2 values (see Table |3, indicating 
that the larger curvature values measured are not statistically 
significant. 

We also verified the possibility to have an additional in- 
trinsic absorption either in the nuclear environment of PKS 




0.5 1 2 

Energy (keV) 



Fig. 2. Best fist spectrum and residuals for the log-parabolic 
model of 0.3-10 keV XRT data (2005, May 22). Data points 
between 0.4 and 0.6 keV were removed to avoid residual cali- 
bration uncertainties in this energy range. 

0548-322 or in the host elliptical galaxy. We recall that the 
exclusion of data points between 0.4 and 0.6 keV do not allow 
an accurate estimation of the absorption. A simple power law 
best fit with a free A^h for the May 22 observation resulted in 
a slightly larger absorption column density of (4.2± 0.6)xl0 20 
cm" 2 . However, the x 1 was found comparable to the fit with 
N H fixed at the Galactic value Ot- 2 /dof= 1.14/408) and, again, 
the residuals showed a clear spectral curvature. The May 22 
spectrum was also fit with the log-parabolic model leaving Nu 
free to vary. The fit was good Ct 2 /dof= 1.03/407) and statis- 
tically equivalent to the case with fixed Galactic absorption. 
We obtained an absorption column density of (2.1+ 1.0)xl0 20 
cm" 2 , in agreement with the Galactic value, and a curvature 
b = 0.21+0.10. 

Therefore, we conclude that the X-ray spectrum of PKS 
0548-322 is intrinsically curved and that it is well described 
by a log-parabolic model with no excess absorption required, 
like the two well known HBL sources Mkn 421 and Mkn 501. 

4. BeppoSAX observations 

In this section we present a re-analysis of BeppoSAX data 
considering the log-parabolic model, not used by 
|Costamante et al. (200 1)| to evaluate the spectral curva- 
ture and compare it to that seen by the XRT in different 
brightness states. 

In the two BeppoSAX observations in 1999 (February 20 
and on April 7) PKS 0548-322 was too faint to be detected 
by the PDS instrument and therefore we analyzed only LECS 
and MECS data. Events for spectral analysis were selected in 
circular regions of 6' and 4' and in the energy bands 0.1-2.0 
keV and 2.0-10.0 keV for the LECS and MECS, respectively. 
Background spectra were taken from the blank field archive at 
the ASI Science Data Center. In Table |2] we reported also the 
log-parabolic best fit parameters of these BeppoSAX observa- 
tions. 

We first note that in 1999 the typical 2-10 keV flux of PKS 
0548-322 was about a factor of 2 lower than in 2005. This 
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Fig. 3. Optical, UV and X-ray Spectral Energy Distribution of PKS 0548-322 from BeppoSAX (1999 February 20, open circles; 
1999 April 07: open squares), Swift XRT (2005 May 22, filled circles) and Swift UVOT (2005 May 21: filled squares) observa- 
tions. The thick dashed lines represent the log-parabolic best fit models to the X-ray data points. The thin dashed line is a power 
law extrapolation of the UV spectral data to the soft X-ray band. Error bars indicate statistical uncertainties at the lcr level. The 
host galaxy contribution has not been subtracted from UVOT data (see Sect. 12. 2\ . 



change of luminosity was accompanied by a different spectral 
curvature (b 0.45) while the peak energy appears more sta- 
ble. To verify that this change of curvature is real we fit the 
BeppoSAX. data with b kept frozen to the value measured by 
Swift in the May 22 observation. We found an increase of x 2 
corresponding to an F-test probability of 4 x 10~ 6 and 2 x 10~ 5 
for the observation of February 20 and April 7, respectively. 
These results make us confident that the change of curvature 
is very significant. We also verified the possibility to have an 
extra-absorption in the local frame of PKS 0548-322. The 
two BeppoSAX 1999 spectra were fit with the log -parabolic 
model with free Nh- In both pointings we found an absorption 
column density of Nh = (3.2 + 1.2) x 10 20 cirT 2 , consistent 
with the Galactic value, and curvatures b = 0.39 + 0.18 and 
b = 0.32+0.20 for the February 20 and April 7 observations, re- 
spectively. Moreover, the fits were statistically equivalent with 
the ones with fixed absorption. We conclude that the X-ray 
spectrum of PKS 0548-322 was intrinsically curved also in 
the two BeppoSAX 1999 observations and is well described by 
a log-parabolic model with Galactic absorption. 



5. Discussion 

Our spectral analysis of the recent series of Swift XRT obser- 
vations of PKS 0548-322 has shown that the X-ray spectrum 
of this BL Lac object is characterized by a significant cur- 
vature that is well fitted by a log-parabolic law with a cur- 
vature parameter and peak energy similar to those of other 



HBL sources (e.g. Massaro et al. 2004a Gio mmi et al. 2 002 
ITramacere et al. 20 06 1. Blustin et al. (2004) summarized some 
historical X-ray data of PKS 0548-322 since its first X-ray 
observation of March 1979 and found that the 2-10 keV flux in 

this source was always within the range (~1.5 4.5)xl0~ n 

erg cirT 2 s . The Swift XRT observations in the spring 2005 
described in this paper confirm this rather stable behavior with 
observed fluxes close to the historical maximum. This relative 
stability at X-ray frequencies of PKS 0548-322 contrasts with 
the behavior of the Mkn 421 and Mkn 501 dFossati et al. 2000; 
Massaro et al. 2004a,b) whose X-ray fluxes varied well over 
an order of magnitude. Of course, the sampling of this historic 
X-ray data set is rather poor and we cannot exclude the oc- 
currence of large outbursts on a time scale of a few months, 
possibly associated with spectral changes. 

Figure shows the SEDs of PKS 0548-322 in the UV 
to X-ray range for the long May 22 pointing and the two 
BeppoSAX observations. The main differences between the re- 
cent high state with respect to the previous data is the well ap- 
parent decrease of the curvature whereas the peak energy is 
around the same values. Moreover, the UV data can be con- 
nected to the low energy X-ray points by a power law inter- 
polation (thin dashed line) suggesting that they come from the 
same component. Simultaneous UV data are not available for 
the BeppoSAX observations: however, a low-energy extrapola- 
tion of the X-ray SEDs indicates that at those epochs the UV 
emission of PKS 0548-322 should have been lower than that 
found in 2005. A recent extensive study of a large set of X-ray 
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observations of Mkn 421 (Massaro et al. 2004a, Tramacere et 
al. 2006b) has shown that this HBL source shows a positive 
correlation between the energy and the flux at the SED peak. 
In the case of PKS 0548-322 we do not observe a significant 
increase of the peak energy with flux. This finding can be an 
indication that statistical acceleration works with different effi- 
ciencies in these two sources. 

PKS 0548-322 is a good candidate for a de- 
tection in the TeV range, as already noticed by 
|Costamante & Ghisellini (2002)| In a Synchrotron Self- 
Compton scenario a brightening at X-ray frequencies often 
corresponds to an enhanced TeV luminosity. It is therefore very 
useful to organize a monitoring program of the X-ray flux in 
the next months to trigger a possible TeV detection. The study 
of the evolution of both synchrotron and inverse Compton 
components in the SED will be very useful to understand the 
physical conditions in the nuclear region and, when compared 
with other TeV BL Lacs, to derive a more general picture 
of this class of sources. Thanks to the wide field of view of 
the BAT instrument, and to its fast pointing capability with 
the UVOT and XRT narrow field telescopes, Swift is the best 
suited satellite to perform such a program. 
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